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Abstract. We have searched for 14.4 keV solar axions or more general axion-like 
particles (ALPs), that may be emitted in the Ml nuclear transition of ^^Fe, by using 
the axion-to-photon conversion in the CERN Axion Solar Telescope (CAST) with 
evacuated magnet bores (Phase I). From the absence of excess of the monoenergetic 
X-rays when the magnet was pointing to the Sun, we set model-independent constraints 
on the coupling constants of pseudoscalar particles that couple to two photons and to 
a nucleon g^^ \ - 1.19 + gl^\ < 1.36 x IQ-^^ GeV-^ for Wa < 0.03 eV at the 95% 
confidence level. 
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1. Introduction 

Quantum chromodynamics (QCD), one of the most profound theories in modern physics 
and nowadays universally believed to be the theory of strong interactions, has one serious 
blemish: the "strong CP problem" . It arises from the fact that the QCD Lagrangian has 
a non-perturbative term (the so-called "6-term") which explicitly violates CP invariance 
in strong interactions. A very credible, and perhaps the most elegant solution to the 
strong CP problem was proposed by Peccei and Quinn in 1977 [1, 2]. It is based on the 
hypothesis that the Standard Model has an additional global U{1) chiral symmetry, now 
known as PQ (Peccei-Quinn) symmetry f/(l)pQ, which is spontaneously broken at some 
large energy scale /a- An inevitable consequence of the Peccei-Quinn solution is the 
existence of a new neutral pseudoscalar particle, named the axion, which is the Nambu- 
Goldstone boson of the broken f/(l)pQ symmetry [3,4]. Due to the U{1)pq symmetry 
not being exact at the quantum level, as a result of a chiral anomaly, the axion is not 
massless and is, more precisely, a pseudo Nambu-Goldstone boson. 

The phenomenological properties of the axion are mainly determined by the scale 
/a and closely related to those of the neutral pion. The axion mass is given by the 
relation 



where z = 0.56 is assumed for the mass ratio of the up and down quarks, while 
~ 92 MeV and = 135 MeV are the pion decay constant and mass respectively. 
Furthermore, the effective axion couplings to ordinary particles (photons, nucleons, and 
electrons) are inversely proportional to /a as well, but they also include significant 
uncertainties originating from some model-dependent numerical parameters. It was 
originally thought that the energy scale of the U{1)pq symmetry breaking is equal to 
the electroweak scale, i.e. /a = few, with /cw — 250 GeV. The existence of the axion 
corresponding to such a scale, known as the "standard" axion, was soon excluded by a 
number of experiments. Despite the failure of the standard axion model, it was possible 
to retain the Peccei-Quinn idea by introducing new axion models which decouple the 
t/(l)pQ scale from the electroweak scale, assuming that /a has an arbitrary value much 
greater than 250 GeV, extending up to the Planck scale of 10^^ GeV. Since the axion 
mass and its coupling constants with matter and radiation all scale as l//a, the axion in 
the models with /a ^ few is very light, very weakly coupled, and very long-lived, which 
makes it extremely hard to detect directly. This is why such models are generically 
referred to as "invisible" axion models and they are still viable. Two classes of invisible 
axion models are often discussed in the literature: KSVZ (Kim, Shifman, Vainshtein, and 
Zakharov) or hadronic axions [5, 6] and DFSZ (Dine, Fischler, Srednicki, and Zhitnitskii) 
or GUT axions [7,8]. The main difference between the KSVZ and DFSZ- type axions is 
that the former do not couple to ordinary quarks and leptons at the tree level. However, 
due to the axion's interaction with photons, there is a radiatively induced coupling to 
electrons present for this type of axions, which is a process of higher order and hence 




(1) 
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extremely weak. As far as the interaction with photons is concerned, once the scale 
/a is fixed, the axion-photon coupling constant Qg^^ for DFSZ axions is also fixed, while 
hadronic axion models suggest different values for g^.^. Consequently, this coupling can 
be either suppressed or enhanced. 

Depending upon the assumed value for /a, the existence of axions could have 
interesting consequences in astrophysics and cosmology. The emission of axions 
produced in the stellar plasma via processes based on their couplings to photons, 
electrons, and nucleons would provide a novel energy-loss mechanism for stars. This 
could accelerate the evolutionary process of stars, and thereby shorten their lifetimes. 
Axions may also exist as primordial cosmic relics copiously produced in the very early 
Universe, which makes them interesting candidates for the non-baryonic dark matter. 

So far the axion has remained elusive after over 30 years of intensive research, 
and none of the direct laboratory searches has been able to yield a positive signature 
for the axion. However, data from numerous laboratory experiments and astrophysical 
observations, together with the cosmological requirement that the contribution to the 
mass density of the Universe by relic axions does not overdose the Universe, restrict 
the allowed values of axion mass to a rather narrow range of 10~^eV < ma < 10"^ eV, 
but with uncertainties on either side. Thus the question of whether axions really do 
exist or the Peccei-Quinn mechanism is not realized in Nature still remains open, and 
the exhaustive search for axions continues. Detailed and updated reviews of the axion 
theory and experiments can be found in [9-11]. 

It is expected that pseudoscalar particles like axions should be copiously produced 
in stars by nuclear reactions and thermal processes in the stellar interior. A powerful 
source of axions would be the Sun in particular. As the closest and the best known 
astrophysical object, it is the source of choice for axion searches. Axions or similar 
axion-like particles (ALPs) that couple to two photons could be produced in the core of 
the Sun via Primakoff conversion of thermal photons in the electric and magnetic fields 
of the solar plasma. Such axions would have a continuous energy spectrum peaked near 
the mean energy of 4.2 keV and dying off above ~10 keV. Most of the experiments that 
have been designed to search for these axions are based on the coherent axion-to-photon 
reconversion in a laboratory transverse magnetic field (the axion helioscope method [11- 
21]), or in the intense Coulomb field of nuclei in a crystal lattice of the detector (the 
Bragg scattering technique [22-26]). 

Due to the axion-nucleon coupling, there is an additional component of solar axions. 
If some nuclei in the Sun are excited either thermally (e.g. those with low-lying levels like 
^^Fe and ^^Kr) or as a result of nuclear reactions (e.g. ''Li* nuclei produced by the ^Be 
electron capture), axion emission during their nuclear de-excitations could be possible. 
Such axions would be monoenergetic since their energy corresponds to the energy of 
the particular nuclear transition which produced them. With axions being pseudoscalar 
particles, the allowed values of angular momentum and parity carried away by them in 
such nuclear transitions are 0~, 1"^, 2~, . . . , which means that axions can be emitted in 
magnetic nuclear transitions. Up to now, these monoenergetic solar axions have mostly 
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been searched for by using the resonant axion absorption process in targets made of 
the same nuchdes in a terrestrial laboratory [27-33], or via processes based on axion- 
electron interactions, like Compton conversion of axion to photon and the axioelectric 
effect [34,35]. 

In this paper we present the results of our search for monoenergetic 14.4 keV solar 
hadronic axions by using the CERN Axion Solar Telescope (CAST) setup. We assumed 
that such particles could be emitted from the Sun by de-excitation of thermally excited 
^^Fe nuclei [36]. This stable isotope of iron (natural abundance 2.2 %) is expected to be 
a suitable emitter of solar axions. It is exceptionally abundant among heavy elements 
in the Sun (solar abundance by mass fraction 2.8 x 10~^). Also, its first excited nuclear 
state {E* = 14.4 keV) is low enough to be relatively easily thermally excited in the hot 
interior of the Sun {kT ~ 1.3 keV). The excited ^^Fe nucleus relaxes to its ground state 
mainly through the emission of the 14.4 keV photon or an internal conversion electron. 
Since this de-excitation occurs dominantly via an Ml transition (E2/M1 mixing ratio is 
0.002), an axion could also be emitted. 

In our attempt to detect ^^Fe solar axions we relied on the axion helioscope method 
by using the latest and currently the most sensitive solar magnetic telescope, CAST, 
which is located at CERN. When its 9.26 m long LHC dipole test magnet is oriented 
towards the Sun, solar axions could convert to photons of the same energy via inverse 
Primakoff process while traversing the 9 T magnetic field produced in the two parallel 
bores inside the magnet. We have searched for a 14.4 keV peak in a spectrum recorded 
by the X-ray detector placed at the far end of the magnet facing the apertures of the 
bores. If observed, this peak could be interpreted as the result of the conversion of 
the ^''Fe solar axions into photons inside the magnet bores, and hence as the direct 
signal for such axions. As a contrast to the previous searches for these axions, which 
relied solely on the axion-nucleon coupling, our search involved not only the axion- 
nucleon interaction (in the emission process) but also the axion-photon interaction (in 
the detection process). This allowed us to explore the relation between axion-photon 
and axion-nucleon coupling constants for a wide range of axion masses. Although we 
focus on the axion because of its theoretical motivation, our results also apply to similar 
pseudoscalar particles that couple to two photons and can be emitted in the nuclear 
magnetic transition. 

2. Axion emission from ^^Fe nuclei in the Sun 

The axion-nucleon coupling arises from two contributions : the tree-level coupling of the 
axion to up- and down-quarks, and a contribution due to the generic axion-pion mixing, 
a phenomenon which is the result of the axion-gluon coupling, and the fact that axion 
and pion are bosons with the same quantum numbers, so they mix. This means that 
even a hadronic axion has a couphng to nucleon, although it does not couple directly to 
ordinary quarks. 
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The effective Lagrangian for the axion-nucleon interaction can be written as 

£aN = ^ a V'n75 (^aN + ^^N^s) V^N , (2) 

where a is the axion field, ~( n ) nucleon doublet, and is the Pauli 

matrix. The isoscalar g^^ and isovector g^^^ axion-nucleon coupling constants are model 
dependent, i.e., they depend on the details of the theory implementing the Peccei-Quinn 
mechanism. For example, in hadronic axion models they are related to the scale /a by 
expressions [37, 38] 

g% = -'^li2S + 3F-D) (3) 

/a 6 

and 

o m-K 1 , ^1 — z 

Here, ttin is the nucleon mass, while the constants F = 0.462 and D = 0.808 [39] are the 
antisymmetric and symmetric reduced matrix elements for the SU (3) octet axial- vector 
currents. They are determined from the hyperon semileptonic decays and flavor SU (3) 
symmetry. The fiavor-singlet axial-vector matrix element 5* is still a poorly constrained 
parameter. It can be estimated from the polarized nucleon structure functions data, 
but suffers from large uncertainties and ambiguity. 

We focused our attention on the decay of the 14.4 keV first excited state of ^^Fe 
nucleus to the ground state via axion emission, a process that competes with ordinary 
Ml and E2 gamma decay. In general, the axion-to-photon emission rate ratio for the 
Ml nuclear transition calculated in the long-wavelength limit is [40] 



r 



a 




\ A;^ / 27ra 1 + 5^ 



^aN/5 + ^^N 



1 2 



(/io-l/2)/3 + /i3-r/ ' 

where fca is the momentum of the outgoing axion, k^ represents the photon momentum, 
and a ~ 1/137 is the fine structure constant. The quantities /io=0.88 and /i3=4.71 
are the isoscalar and isovector nuclear magnetic moments respectively, given in nuclear 
magnetons. The parameter 5 denotes the E2/M1 mixing ratio for the particular nuclear 
transition, while /5 and rj are nuclear structure dependent ratios. Their values for 
the 14.4 keV de-excitation process of an ^^Fe nucleus are 5=0.002, f] = —1.19, and 
rj = 0.8 [36]. Using these values in equation (5) we find 

^ = l.82i-1.19gl^ + gl^y. (6) 

In the above expression we made approximation (k^/k^)^ ~ 1 for the phase space factor 
ratio. This is based on the assumption that the axion mass is negligible compared to 
the axion energy which equals, in our case, the transition energy of 14.4 keV. Since the 
coupling constants g% and (/aN model dependent, we can consider the parameter 
g"^ = (— 1.195'aN + daTSs) ^ uukuowu parameter characterizing not only the 
axion-nucleon coupling but also, more generally, the nucleon coupling to any axion- 
like particles that could be emitted in the Ml nuclear transition. In terms of the axion 
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mass, equation (6) can be expressed by combining equations (1), (3) and (4) as 

I^^3.96xl0-(^)\ (7) 

where S = 0.4 [41] is assumed. 

To estimate the flux of monoenergetic solar axions emitted from ^^Fe nuclear 
de-excitations, the calculation was performed as in [36,42]. The ^^Fe nucleus can 
be thermally excited to its first excited state in the Sun since its excitation energy, 
i?*=14.4 keV, is comparable to the Sun's core temperature {kT ~ 1.3 keV). The 
probability for this thermal excitation is given by the Boltzmann distribution 

(2 Ji + 1) e-^y^^ 

(2Jo + l) + (2Ji + l)e-^*/'=^' ^ ' 

where Ji = 3/2 and Jq = 1/2 are total angular momenta of the first excited and 
ground states respectively. In further calculations we can neglect the second term in 
the denominator because q-^*/^'^ <^ i. The 14.4 keV axion emission rate per unit mass 
due to the de-excitation of thermally excited ^''Fe nuclei in the Sun is 

A4=Ar 1^,11^ g-^s-^ (9) 

where A/" = 3.0 x 10^^ g~^ is the number of ^^Fe nuclei per 1 g of solar matter [43], and 
= 1.3 X 10"^ s is the mean lifetime of the first excited nuclear state of ^^Fe associated 
with the partial gamma decay width of that state. Owing to Doppler broadening 
caused by the thermal motion of ^^Fe nuclei in the hot solar interior, the axion emission 
spectrum is a Gaussian with the standard deviation parameter 




a{T) = E*d—, (10) 



where T denotes the temperature at the location in the Sun where the axion is produced, 
while m is the mass of the ^^Fe nucleus. This spectrum is approximately centered at 
the transition energy Ei = E* = 14.4 keV because the axion energy shift, caused 
by the recoil of the ^^Fe nucleus in the emission process (~ 1.9 x 10~^ eV) and the 
gravitational redshift due to the Sun (~ 1.5 x 10~^ eV), is negligible compared with 
the Doppler broadening of the spectrum (FWHM= 2.35 cr(T) ~ 5 eV). Following these 
arguments, we can write the differential ^'^Fe solar axion flux expected at the Earth as 

d^.jE,) ^ J_ ^ 1 r (^a-^i) 

dE, Audi Jo " v^a(r(r)) ''''^ [ 2a(T(r))^ 

X p(r) 47rr^ dr . (11) 

Here E^^ is the axion energy, ^0 denotes the average distance between the Sun and the 
Earth, Rq is the solar radius, while T(r) and p(r) are the temperature and the mass 
density in a spherical shell with the radius r in the solar interior, respectively. 

The total flux of ^^Fe solar axions at the Earth can be calculated by integrating the 
differential axion flux (equation (11)) with respect to the axion energy. This leads to 



4. . I'^iE. 

dE„, 
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with the help of equations (8) and (9). As in the framework of the standard solar model, 
we assumed that the number of ^^Fe nuclei per unit mass in the Sun is uniform, i.e., 
that jV is independent of r. Using the standard solar model data for the temperature 
and mass density distributions in the Sun as functions of the radius r [44] , we evaluated 
the integral in the above expression and found that 

<l>a = 4.56 X 10"=^ {g'^y cm-2 . (13) 

The corresponding solar axion luminosity is calculated to be | 

= 7.68 X 109(^7:^)2 Lq, (14) 

where Lq = 3.84 x 10^^ W is the solar photon luminosity. Arguments related to the 
measured solar neutrino flux (La < 0.1 Lq) [46] lead to 

gll\< 3.6x10-'. (15) 

3. The detection of solar axions in the CAST experiment 

The CERN Axion Solar Telescope (CAST) experiment is the most recent 
implementation of the axion helioscope technique [12]. It searches for solar axions 
and similar particles with unprecedented sensitivity in the sub-eV mass range which is 
comparable to the astrophysical constraints on these particles. Like most of the axion 
search experiments in the past 30 years, this experiment relies on the axion coupling to 
two photons, a generic property of axions and ALPs given by the effective Lagrangian 

1 

C^^ = --g^^F'"'F^^a = g^^E-Ba, (16) 

where a is the axion field, F'^'^ the electromagnetic field strength tensor, and F^^ its 
dual. This interaction can also be expressed in terms of electric E and magnetic B field 
of the coupling photons, as shown in the above expression. The effective axion-photon 
coupling constant gg^y is given by 



a 



E_ 2 (4 + ^) 
iV ~ 3(l + z) 



a f E , , , 

— - 1.95 ±0.08 . (17) 



27r/a \N 

Here E and N are the model dependent coefficients of the electromagnetic and color 
anomaly of the axial current associated with the f/(l)pQ symmetry, respectively. In 
DFSZ axion models their ratio is fixed to E/N = 8/3, while for the hadronic axions this 
ratio can take different values depending on the details of each model. 

Axions and similar particles with a two-photon interaction of the form given by 
equation (16) can transform into photons, and vice versa, in external electric or magnetic 
fields. Therefore, these particles could be produced in stars by the Primakoff conversion 

I Following the calculations in [45] wc estimated that processes like the Compton process on nuclei 
and proton-Helium bremsstrahlung, that arc also based on the axion-nucleon coupling, are suppressed 
relative to the hadronic axion emission from the Ml nuclear transition of ^^Fe. 
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of thermal photons in the Coulomb fields of nuclei and electrons in the stellar plasma. 
On the other hand, the reverse process has served as the basis for various experimental 
methods to search for these particles. The axion helioscope technique [12] is one such 
method. The essence of this idea is to search for solar axions using a long dipole 
magnet in a laboratory. Inside the magnet, while it is oriented towards the Sun, the 
incoming axion couples to a virtual photon provided by the transverse magnetic field 
and converts into a real photon via the reverse Primakoff process a + 7virtuai 7- This 
photon has the energy equal to the axion energy and can be detected with a suitable 
X-ray detector placed at the far end of the magnet opposite the Sun. Assuming that 
there is a vacuum in the conversion volume (i.e. magnetic field region), the axion-photon 
conversion probability is [14] § 

I 2 ) q^L^ ' I 2 
where B is the magnetic field, L is the length of the conversion volume in the direction of 
the incoming axion propagation, and q = m1/2Ea represents the momentum difference 
between the axion and the photon of energy E^. The conversion probability is maximal 
if the axion and photon fields remain in phase over the length of the conversion region, 
i.e., when the coherence condition gL < vr is satisfied [14]. This restricts the sensitivity 
of a helioscope to a specific range of axion masses, e.g., for a 10 m long magnet and 
axion energy of ~10 keV the coherence condition sets the limit of ma ^0.03 eV on the 
axion mass, up to which such an experiment is sensitive. However, coherence can be 
maintained for higher axion masses if the conversion volume is filled with a buffer gas 
such as hehum [13]. In this case, photons acquire an effective mass m-y whose value is 
determined by the gas pressure, and the axion-photon momentum difference becomes 
q = — /2Es,. As a result, the coherence is restored for a narrow axion mass 
range ^ m^, where the effective photon mass matches the axion mass. 

The first experiment to use the axion helioscope technique was performed at BNL in 
the early 1990s [14]. Following this experiment, the Tokyo Axion Helioscope continued 
the search for axions using the same method with much improved sensitivity [15-17]. At 
present, the most sensitive axion helioscope is the CERN Axion Solar Telescope (CAST). 
The main component of CAST is a decommissioned prototype of a twin aperture LHC 
superconducting dipole magnet, which serves as a magnetic telescope. It provides a 
transverse magnetic field of 9.0 T inside the two parallel, straight, 9.26 m long bores. 
The aperture of each bore is 43 mm, so the total cross-sectional area is 2 x 14.5 cm^. In 
terms of the parameter [BL)"^, which according to equation (18) determines the axion- 
photon conversion probability, the CAST magnet is ~ 80 times more efficient as an 
axion-to-photon converter than the best competing helioscope in Tokyo. To optimize 
cast's performance, the magnet is installed on a moving platform which allows it to 
track the Sun ±8° vertically and ±40° horizontally. Thus, it can be aligned with the Sun 
for approximately 1.5 h both during sunrise and sunset every day throughout the year. In 

S We use natural units with h = c = 1. 
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order to detect photons coming from the magnet bores, as a result of axion conversion in 
the magnetic field, several low-background X-ray detectors are installed on both ends of 
the magnet. Until 2007, a conventional Time Projection Chamber (TPC) was located 
at one end, covering both magnet bores, to detect photons originating from axions 
during the tracking of the Sun at sunset. It was then replaced by two MICROMEGAS 
detectors, each attached to one bore. On the other side of the magnet, there is another 
MICROMEGAS detector covering one bore, and an X-ray mirror telescope with a pn- 
CCD chip as the focal plane detector at the other bore, both intended to detect photons 
produced from axions during the sunrise solar tracking. More details about the CAST 
experiment and detectors can be found in [18-21,47-49]. 

To cover as wide as possible range of potential axion masses, the operation of the 
CAST experiment is divided into two phases. During the Phase I (2003-2004) [19,20] 
the experiment operated with vacuum inside the magnet bores and the sensitivity was 
essentially limited to ma < 0.02 eV due to the coherence condition. In the second phase 
(so-called Phase II) which started in 2005, the magnet bores are filled with a buffer gas 
in order to extend the sensitivity to higher axion masses. In the first part of this phase 
(2005-2006) ^He was used as a buffer gas. By increasing the gas pressure in appropriate 
steps, axion masses up to ~0.4 eV were scanned and the results of these measurements 
supersede all previous experimental limits on the axion-photon coupling constant in this 
mass range [21]. To explore axion masses above 0.4 eV, ^He has to be used because it 
has a higher vapor pressure than ^He. This allows us to further increase gas pressure 
in the magnet bores and to reach axion masses up to about 1 eV in the ongoing second 
part of Phase II that started in 2007 and is planned to finish by the end of 2010. 



4. Measurement and data analysis 

The CAST experiment is primarily designed to search for axions or axion-like particles 
that could be produced in the Sun by the Primakoff conversion of thermal photons. 
Their expected energy spectrum [20], 

^^liEA = 6.02 X 10^° gl El-''' e-^^l^-^^^ cm-' s'' keV'^ , (19) 
dE„ 



■'a 



(where energies are in keV) has the peak at 3 keV, mean energy of 4.2 keV, and vanishes 
above 10 keV. Since the conversion of these particles inside the CAST magnet bores 
would produce photons of the same energies, the X-ray detectors used in CAST are 
optimized for the efficient detection of photons in the 1-10 keV range. To search for 
14.4 keV photons that might originate from the conversion of the ^'^Fe solar axions we 
used only the data provided by the TPC detector because the other detectors have very 
low sensitivity to photons with energies above ~10 keV. 

The CAST TPC incorporates the well-known concepts of drift chambers and Multi- 
Wire Proportional Chambers (MWPC). It has an active volume of 10 x 15 x 30 cm^ filled 
with an Ar (95%) -|- CH4 (5%) gas mixture at atmospheric pressure, where the incoming 
particle interacts with the gas producing free electrons. The detector's volume is 10 cm 
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long in the drift direction that is parallel to the magnet axis, while the 15 x 30 cm^ 
cross section, covering both magnet bores, is perpendicular to this direction. The drift 
region is bounded on the front side with a 15 x 30 cm^ drift electrode, biased at -7 kV, 
that is made of a thin aluminum layer and covers the entire inner side of the chamber 
wall closest to the magnet. On the opposite end, i.e., back side of the chamber, there 
are 3 planes of sense wires parallel to the drift electrode: one anode plane at +1.8 kV 
with 48 wires placed between two grounded cathode planes with 96 wires each. The 
wires in both cathode planes are perpendicular to the anode wires. The spacing between 
two adjacent wires of the same plane is 3 mm. The gap between the anode and the 
inner cathode plane (the one closest to the drift electrode) is 3 mm, while the distance 
between the anode and the outer cathode plane is 6 mm. This asymmetric configuration 
enhances the induced signals on the cathode wires in the inner plane, which are the ones 
being read out by the front-end electronics, together with the anode wires. Each of these 
wires is read out individually, and since the anode wires are perpendicular to the cathode 
wires, a very good two-dimensional position resolution can be obtained. This allows us 
to distinguish spatially localized photon events that may represent the axion signal from 
the background events characterized by the long tracks of cosmic rays. The detector 
is connected to the magnet by two 6 cm diameter entrance windows made of a thin 
aluminized mylar foil stretched on a metallic grid. These windows allow the photons 
coming from the magnet bores to enter the detector. For a detailed description of the 
TPC detector, its shielding, front-end electronics, and data acquisition system, we refer 
to [49,50]. 

To characterize the CAST TPC detector, a series of calibration measurements was 
performed using the X-ray beams with very accurately calibrated energies and intensities 
at the PANTER X-ray facility of MPE in Munich [49]. Photons originating from the 
axion conversion in the CAST magnet would enter the detector only through its two 
entrance windows, their direction being parallel to the magnet axis, i.e., perpendicular to 
the anode and cathode wire planes. Such a photon would normally produce a point-like 
energy deposition in the region of the detector's sensitive volume facing the windows. 
From the analysis of the PANTER data it was realized that the pulses induced on the 
sense wires (hits) in such event are clustered on several contiguous anode wires (so- 
called anode signal cluster), as well as on several contiguous cathode wires (cathode 
cluster). This characteristic profile of a photon event can be seen in the example shown 
in figure 1. The left and the right plot display the time evolution of the pulse induced 
in each anode and each cathode wire, respectively, as recorded by the flash-ADCs. 
By studying the cluster properties and topology, and exploiting the two-dimensional 
position reconstruction capability of the detector, we were able to establish a set of 
rules (cuts) to distinguish genuine photon events that might be an axion signal from 
background events coming from cosmic rays and natural radioactivity. These cuts take 
into consideration the total number of anode and cathode clusters in the event, their 
time correlation, cluster multiplicity given by the number of hits in each cluster, and the 
position of clusters in the anode and cathode wire planes with respect to the entrance 
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Figure 1. Time evolution of the pulse induced in each wire as recorded by the flash- 
ADCs in a photon-like event. Left (right) plot shows pulses in anode (cathode) wires. 

windows of the detector [49]. We should emphasize that these cuts were derived from 
the data obtained during cahbration measurements at PANTER where monochromatic 
photon beams with energies from 0.3 to 8 keV were used. This energy range was chosen 
because it covers the bulk of the expected spectrum of solar axions produced by the 
Primakoff process, and these axions are, as stated before, the main subject of CAST's 
search for solar axions. It is expected that the similar cuts can be applied to single out 
14.4 keV photon events that could indicate the conversion of ^^Fe solar axions inside 
the bores of the CAST magnet. In order to verify this assumption, we made additional 
calibration measurements with a ^^Co source placed in front of each of the two TPC's 
windows. ^'^Co was the source of choice for this task because in the course of the ^^Co 
decay to the ground state of ^^Fe, several photons are emitted, and one of them has the 
energy of 14.4 keV. 

Figure 2 shows the spectrum of ^''Co source after the selection cuts were applied 
to the TPC data. As a starting point to reject events that were not produced by 
photon interactions in the detector's volume, we followed the information obtained from 
the PANTER calibration measurements and set the requirement that only events with 
one anode cluster and one cathode cluster should be considered as relevant photon 
events. This condition also allowed us to straightforwardly match these two clusters in 
order to obtain the two-dimensional position of the point- like event. Our study of the 
cluster properties in the events that constitute the 14.4 keV peak in the ^^Co calibration 
spectrum resulted in a set of additional cuts given in table 1. The first two cuts consider 
the spread of the clusters. The spread of the anode cluster is due to the diffusion of 
the electron cloud along the drift distance from the interaction point of an incoming 
photon to the anode wires plane. The larger spread of the cathode cluster is due to the 
development of the avalanche process along the anode wires. The "time correlation" cut 
reflects the fact that in a real photon event both anode and cathode cluster originate 
from the same avalanche process induced in the proximity of the anode wires by the 
initial ionization electron cloud. The purpose of this cut is to reject events with spurious 
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Table 1. List of software cuts applied to the CAST TPC data in the search for ^^Fe 
solar axions. 

Cut Condition 



Anode cluster multiplicity 2 or 3 hits in the anode cluster 

Cathode cluster multiplicity 2 to 8 hits in the cathode cluster 

Anode-cathode cluster time correlation Time difference between anode and cathode cluster 

must be in the range —0.15 to 0.02 /is 
No saturation No hit in the cathode cluster reaching the upper limit 

of the flash- ADC dynamical range 
Position 2-D coordinates of the event should be inside 

the area facing the magnet bores 



03 




ADC units 

Figure 2. TPC calibration spectrum measured with '^^Co source after the off-line 
analysis cuts were applied to the data. Two peaks arc present: Fc X-rays (Kq,) at 
6.4 keV and 14.4 kcV photons corresponding to the 7-transition in '^^Fe. The grey line 
shows the combined fit of the peaks plus background while arrows indicate the region 
of interest for the axion signal. 



clusters that might be produced by the effects of noise in contiguous wires and mimic the 
real photon events. The "no saturation" cut is related to the fact that energy deposited 
in each event is calculated using the total strength of the cathode cluster obtained by 
adding up the pulse height of every hit in the cluster. If any of the hits has a pulse 
higher than the one that the fiash-ADC can handle, the calculated energy would be 
incorrect, and thus such event should be rejected. Finally, the "position" cut is used to 
reject events with two-dimensional positions outside the region where the axion signal is 
expected, i.e., out of the area facing the two magnet bores. The application of all these 
cuts in the off-line analysis of data recorded during normal CAST operation reduces the 
background by approximately two orders of magnitude with respect to the raw trigger 
rate without significantly reducing the efficiency of the detector. 
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The hardware efficiency of the TPC to detect 14.4 keV photons was estimated 
by Monte Carlo simulations, using the GEANT4 software package. The geometry of 
the detector, the transparency of the entrance windows, as well as the opacity of the 
gas inside the sensitive detector volume were taken into account, giving the efficiency 
of 13%. Here we note that the validity of the simulations was checked by comparing 
the computed efficiencies with those obtained from the PANTER measurements [49] for 
photon energies of 6.4 and 8 keV (the two highest energies at which measurements were 
performed). The agreement between the corresponding efficiencies is within 4 and 2 
percent respectively. Therefore, we regarded the efficiency computed by simulation as 
a good approximation for the TPC hardware efficiency at 14.4 keV. The additional 
efficiency loss of 18% in the off-line analysis, due to the software cuts applied to 
reduce the background, was calculated using ^^Co calibration data. Hence, the overall 
(hardware + software) TPC detection efficiency of 10.7% for 14.4 keV photons that may 
come from conversion of solar axions was obtained. Finally, data from ^''Co calibration 
runs were also used to determine the energy resolution of the TPC at 14.4 keV which 
defines the energy region of interest in our search for the signal of ^^Fe solar axions. It 
was found that the width of the 14.4 keV peak is characterized by the standard deviation 
parameter a = 1.77 keV, giving a full width at half maximum of 4.16 keV. 

For the study presented in this article, we used the data acquired with the TPC 
over the period of ~5.5 months in 2004 during Phase I of the CAST experiment. The 
energy response of the detector was calibrated periodically using the ^^Fe X-ray source. 
No systematic shift in the energy scale with time was found. Since the calibration runs 
took place every six hours, we were able to characterize, with very good precision, small 
gain variations due to the fluctuations of environmental parameters that affected the 
gain of the detector, and to correct the calculation of the energy deposited in each event 
accordingly. 

Data quality checks were performed both on-line (during data acquisition) and 
off-line by using the "slow control" data that were continuously recorded during the 
operation of the experiment to monitor various experimental parameters (e.g. the 
magnetic field strength, magnet pointing direction with respect to the current position 
of the Sun, magnet temperature, various pressures, etc.). As a result, data qualified 
for the analysis were accumulated during 2819 effective hours, out of which 203 hours 
correspond to the tracking data (data acquired in the "axion-sensitive" conditions, i.e., 
while the magnet was tracking the Sun), and 2616 hours correspond to the background 
data (data acquired during the non-tracking periods). The background data were used to 
estimate and subtract the true background contribution in the tracking data spectrum. 
Therefore, these data were taken while the magnet was parked in well-defined positions 
close to where it was passing during the Sun tracking, but at times when the Sun was 
not in view. 

The energy spectrum of the events reconstructed from the tracking data that passed 
all the software cuts in the off-line analysis, as well as the measured background spectrum 
after applying the same cuts, are shown in figure 3 on the left side. Both spectra are 
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properly normalized so they can be compared. The vertical dashed lines delimit the 
interval where the peak coming from the excess of 14.4 keV photon events, that would 
indicate a signal for ^''Fe solar axions, is expected. As stated earlier, the TPC energy 
resolution (a) at the 14.4 keV peak is 1.77 keV. Therefore, we regarded this 14.4±3.5 keV 
interval as the ±2(j (95.45%) signal region. To extract the signal, we subtracted the 
background spectrum from the tracking one, and the resulting signal spectrum can be 
seen on the right side in figure 3. The number of detected photon events in the signal 
region is —71 ± 57. Consequently, only an upper limit on any axion signal expected in 
this energy range could have been set. 

Since the energy resolution of the TPC detector does not allow us to look at the 
very narrow 14.4 keV line alone, one also picks up a fairly broad part of the Primakoff 
solar axion spectrum which could mask the yield of ^^Fe solar axions. This means that 
our method will not be able to identify their contribution for certain combinations of 
axion-photon and axion-nucleon coupling constant values. To exclude cases where the 
11-18 keV tail of the Primakoff flux exceeds the ^^Fe solar axion flux in the measured 
signal spectrum, we set a constraint on the Primakoff axion emission by requiring 

1-18 keV d^^(E ) 

/ < ■ (20) 

J 11 koV a£/a 

Using equations (13) and (19) this bound translates to 

g^,<lMxlO'''g:^ GeY-\ (21) 

and thus restricts the region of g^,^ — g^ parameter space in which our method, 
complemented by the use of the TPC detector, is strongly sensitive to ^^Fe solar axions, 
as can be seen on the right side in figure 4 (region below the "Det" line). In this view, 
our analysis involves searching for the Gaussian-like energy spectrum of axion-induced 
photons, 



Nf,{E) = — ^iVgexp 



9^2 ' (22) 

that describes the shape of the expected 14.4 keV axion peak in the measured signal 
spectrum. The analysis was performed by standard minimization. The energy 
resolution a = 1.77 keV and the position of the peak E'a = 14.4 keV were fixed during 
the fitting. With t = 203 hours of Sun tracking and A = 2x 14.5 cm^ cross-sectional area 
of both magnet bores, the minimum of x^/d.o.f.=ll. 47/13 corresponds to the number of 
signal events = —42 ± 27. This result is consistent with zero, thus giving no evidence 
for ^''Fe solar axions in our search. Therefore, an upper limit on the number of signal 
events, TVs < 32 counts at 95% C.L., was obtained following the Bayesian approach [9] 
by calculating the value of that encompasses 95% of the physically allowed part 
{Ns > 0) of the Bayesian probability distribution. The corresponding fits are plotted in 
figure 3 (right side). The effect of the Primakoff axion tail on the result of our analysis 
was studied in the parameter space regions where the ^''Fe solar axion contribution is 
comparable to or significantly less than the Primakoff one and discussed in section 5. 
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Figure 3. Left: Energy distribution of events recorded during 203 hours of Sun 
tracking (• ), i.e. while the TPC was sensitive to solar axion signals, compared to 

the background spectrum ( ) measured during non-tracking periods in total of 

2616 hours. Right: The subtracted spectrum (• ) together with the expectation for 

the best fit iVg ( ) and for the 95% C.L. hmit on TVs ( )• The signal region 

(14.4 ± 3.5 keV) that is expected to cover the 14.4 keV peak containing photon events 
coming from the conversion of ^^Fe solar axions in the magnet bores is also shown in 
both plots. 



5. Results and discussion 

The number of photons expected to be detected by the TPC, coming from the magnet 
bores as a result of the conversion of ^''Fe solar axions in the magnetic field, is 

Ns = j P.,{E.) A t e{E,) dE, , (23) 

where e(-E'a) is the TPC detection efficiency for the photons of energy E^,. We can 
rewrite the axion-photon conversion probability P^^{E^, given by equation (18), in a 
more convenient form as 

P.,(E.) = 1.736 X lO^f^Vf— Vf^V^ sin^ fi^) .(24) 
' VGeV~V V9Ty V9.26my \2 J ^ ' 

One should recall that the conversion probability depends implicitly on the axion 

mass ma and the axion energy E^. through the axion-photon momentum difference 

q = m\/2Es,. However, since the ^'^Fe solar axions are nearly monoenergetic, which 

is the result of their energy spectrum being very narrow (FWHM ~ 5 eV), it is a 

reasonable approximation that the conversion probability is constant over their energy 

range and equal to = -Pa7(14.4 keV). Hence, the equation (23) can be written as 

iVs = $a P.-y A t ei4.4 , (25) 

where €14,4 = 0.107 is the detection efficiency of the TPC at 14.4 keV. Combining 
equations (24) and (13), equation (25) in our case becomes 

iVs^ 1.728 X 10- toS)^(^)^^s,„^(«^). (26) 

where for the magnetic field we used its effective value of 8.83 T, due to the fact that 
during 82% of the data taking time the magnetic field was 8.79 T, while for the rest of 
the time its value was 9 T. 
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Figure 4. Left: The upper limit (95% C.L.) on the product Qaj g'^§s as a function of 
the axion mass ma, imposed by the CAST's search for ^^Fe solar axions. Right: The 
upper Hmit on g^-y versus g^^, based on the relations ga-y gl§i < 1-36 x 10"^^ GcV~-^ 
and ga-y < 3.5 x 10~^" GcV~^ for TOa<0.03 eV, is shown. The boundary denoted as 
Lum-Fe is a constraint due to the ^^Fc solar axion luminosity, while the Dot line is 
due to the detector resolution. This line divides the 3a7~5aN parameter space into two 
regions where, roughly speaking, our method is dominantly sensitive to the ^^Fe solar 
axions (below the line) and to the Primakoff axions (above) , while in its proximity the 
sensitivities are comparable. Upper limits from the Primakoff solar axion luminosity 
(Lum-P) [46] and CAST's search for the Primakoff solar axions (CAST-P) [20], that 
rely solely on ga-y, are also displayed for comparisons. 



According to equation (26), the upper limit on the number of signal events A^s < 32, 
resulting from the non-observation of the signal above background in our search for ^^Fe 
solar axions, can be translated into the limit on the product of the parameters Qs,^ and 

cff II 



^'^^^f< 1.36x10-16 



(95%C.L.). (27) 



GeV-i ■ ^4sin2(gL/2) 

We can consider this model-independent limit to be a function of the axion mass 
TTia, as shown on the left side in figure 4. In the mass range ma < 0.03 eV, where the 
axion-photon conversion process is coherent and has maximum probability, the limit is 
mass-independent and its value (95% C.L.) is 

g,,gl^< 1.36x10-'' GeV-' . (28) 

For higher axion masses, an increase of the momentum mismatch q causes a loss of the 
coherence of the axion-to-photon conversion, and suppresses the conversion probability 
due to the factor sin^(x)/x^ in equation (18), with x = qL/2. As a consequence, the 
sensitivity of the experiment to 14.4 keV axions diminishes rapidly with the increase 
of axion mass above ~0.03 eV, thus providing weaker, mass- dependent limit, as can be 
seen from the plot in figure 4 (left side). 

II In the following we always mean j^aNl when we write .gaN- 
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Table 2. Results of analysis in the regions of .9a7^.9aN parameter space where our 
method is sensitive to ^^Fe axions, ^^Fe+Primakoff axions, and Primakoff axions (from 
top to bottom). 



Test 


(5a7 5aN)bcstfit±^ 


(5a7)bcstfit ± a- 




5a7 5^^ (95%C.L.) 


5a7 (95% C.L.) 


hypothesis 


(10-32 GeV-2) 


(10-38 GeV-4) 




(IQ-if^ GeV-i) 


(10-10 GeV-i) 


iVpc 


-2.4± 1.6 




11.47/13 


1.36 






-1.7± 1.7 


-0.5 ± 1.3 


11.44/12 


1.5 


3.8 


iVp 




-1.1 ± 1.1 


12.89/13 




3.5 



The limit from equation (28) allowed us to set the constraint on the axion-photon 
coupling constant as a function of the parameter g'^ for axion masses ma < 0.03 eV as 

1.36 X 10-16 

< ^ GeV . (29) 

fl'aN 

This limit, denoted as "CAST-Fe", is shown on the right side in figure 4. The vertical 
bound labeled "Lum-Fe" at = 3.6 x 10-^ is due to the requirement that the 
axion emission from ^''Fe nuclei in the Sun should not exceed 10% of the solar photon 
luminosity, as was explained in section 2. The method we used to search for ^^Fe solar 
axions can be applied in the region of gi^^- g^ parameter space given by equation (21) 
(i.e. below the line denoted as "Det"), where the ^^Fe solar axion flux exceeds the 
tail of the Primakoff solar axion flux in the energy range of the expected ^''Fe axion 
signal. In this region, the dark brown area is excluded due to the ^^Fe solar axion 
luminosity constraint and the "CAST-Fe" line given by equation (29), while the white 
area is allowed according to our search. Due to the relatively low energy resolution of 
the TPC detector, the CAST search for ^''Fe solar axions is not significantly sensitive 
in the parameter space region (grey) above the "Det" line, where the Primakoff axion 
contribution dominates, while in the proximity of this line we expect that these two 
contributions are comparable. To examine how the Primakoff axion tail affects our 
result (28), we performed two additional analysis of the experimental signal spectrum 
using as fit functions: i) A'pe + A^p, and ii) Ap. Here Ape is given by equations (22) and 
(26), while A^p corresponds to the expected spectrum of the Primakoff axion tail (19) 
multiplied by the axion-photon conversion probability (24) and the detection efficiency. 
The results of these fits are shown in table 2. The suppression of our method's sensitivity 
to the ^^Fe axions in the region above the "Det" line due to the Primakoff axions 
contribution resulted in the upper limit of g^^-y < 3.5 x 10-^° GeV-^ at 95% C.L., which 
is displayed as a red horizontal line in figure 4 (right side). Also, according to table 2, 
the effect of the Primakoff axion tail on our upper limit fi'a7 fl'aN < 1-36 x 10-^^ GeV-^ 
is ~ 10% in the proximity of the "Det" line, while the ^^Fe axions affect the limit 
c/a7 < 3.5 X 10-1° GeV-1 for about 9%. 

The experimental systematic uncertainties on the present limits were studied. 
Regarding the background determination, the null hypothesis test (in areas of the TPC 
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Figure 5. Exclusion plot (95% C.L.) in the g^-y versus ma plane. The limit achieved by 
our research is presented for two values of the parameter reflecting the luminosity 
and detection restrictions (from bottom to top, respectively). The yellow band 
represents typical theoretical axion models with \E/N — 1.95| in the range 0.07-7. 



detector where no signal is expected) was used in order to estimate the systematic 
uncertainty induced by possible uncontrolled dependencies of the background on time, 
position or other experimental conditions. These effects were considered by varying 
the background level until the null hypothesis test yielded a result with a probability 
smaller than 5%. If taken as an uncertainty, this range corresponds to 



+8% 
-10% 



variation 



of the upper limits. As stated earlier, the deviation of our result due to the inclusion of 
possible contribution of the Primakoff axion tail in the fitting procedure was estimated 
to be less than +10%. Other effects, such as uncertainties of the magnet parameters, 
are negligible, while the TPC efficiency affects the upper limit less than -2%. Therefore, 
we estimated, using the quadratic sum of the individual contributions, that the overall 
effect of systematic uncertainties on our upper limit of g^^yg^N ^^^s than . 

It is important to stress that the only axion properties we relied on in the entire 
procedure that led us to equation (27) were that its couplings to photons and nucleons 
have a general form given by the Lagrangians in equations (2) and (16). We did not 
use any specific details regarding the coupling constants g^^-y, g^^ and g^^^ from any of 
the axion models. Therefore, we can consider these coupling constants as free unknown 
parameters that characterize the couplings of axions or general axion-like particles to 
two photons and a nucleon. This allowed us to use equation (27) in order to set the upper 
limit on g^,^ as a function of ma for various values of g'^. Figure 5 shows the exclusion 
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plots of (y'a7 versus ma obtained for two values of the parameter in comparison with 
the "axion models band", i.e., the region of g^^^- m.^ values expected from typical axion 
models with \E/N — 1.95| in the range 0.07-7. The axion-nucleon couplings can vary 
from 3.6 x 10^^ to 3.6 x 10~^ reflecting constraints due to the ^''Fe solar axion luminosity 
and detection sensitivity, respectively. The presented contours do not enter the range 
of parameters, indicated by the yellow band in figure 5, that is predicted by plausible 
axion models to be the best-motivated region to search for axions. Thus, they should 
rather be considered as the limits on the two-photon coupling of axion-like particles that 
are somewhat lighter for a given interaction strength than it is expected for axions. 

These two contours serve as an example to show how our result, given by the 
equation (27), can be used to scale the excluded region in the g^^^/- ma parameter space 
for various choices of g^^ and g^^. In such a manner, it can be generally applied to impose 
the constraints on exotic light pseudoscalar particles that could be emitted in the nuclear 
magnetic transitions and couple to two photons. CAST has also performed a similar 
search for high-energy solar axions and ALPs from ''Li (0.478 MeV) and D(p,7)^He 
(5.5 MeV) nuclear transitions, and the results are reported in [51]. 

6. Conclusion 

The ongoing CAST experiment is primarily designed to search for hadronic axions or 
more general ALPs of continuous energy spectrum, with an average energy of 4.2 keV, 
that could be produced abundantly in the solar core by the Primakoff conversion of 
thermal photons in the electric fields of charged particles in the plasma. Since the 
reconversion of these particles inside the CAST magnet bores would produce photons 
of the same energies, the X-ray detectors used in CAST are optimized for the efficient 
detection of photons in the 1-10 keV range. Here we explored the relation between the 
coupling constants of pseudoscalar particles that couple to a nucleon and to two photons 
by using the CAST setup during the Phase I to look for 14.4 keV monoenergetic solar 
axions and ALPs that may be emitted in the Ml nuclear transition of ^^Fe. The signal 
we searched for, i.e., an excess of 14.4 keV X-rays when the magnet was pointing to the 
Sun was not found, and we set model-independent limits on the coupling constants of 
5'a7 I — 1-19 5?°^ + 91^1 < 1-36 X 10"-^^ GeV~-^ at the 95% confidence level. As a contrast 
to other experiments sensitive on the gl^gl-^ couplings [35,40,52,53] that put some 
constraints in the ~ 10^ — 10^ eV axion mass range, we explored the low mass region 
up to 0.03 eV. 

Acknowledgments 

We thank CERN for hosting the experiment and for the contributions of J. P. Bojon, 
F. Cataneo, R. Campagnolo, G. Cipolla, F. Chiusano, M. Delattre, A. De Rujula, 
F. Formenti, M. Genet, J. N. Joux, A. Lippitsch, L. Musa, R. De Oliveira, A. Onnela, 
J. Pierlot, C. Rosset, H. Thiesen and B. Vullierme. We acknowledge support from 



Search for 14-4 keV solar axions with CAST 



21 



NSERC (Canada), MSES (Croatia) under the grant number 098-0982887-2872, CEA 
(France), BMBF (Germany) under the grant numbers 05 CC2EEA/9 and 05 CClRDl/0, 
the Virtuelles Institut fiir Dunkle Materie und Neutrinos - VIDMAN (Germany), GSRT 
(Greece), RFFR (Russia), the Spanish Ministry of Science and Innovation (MICINN) 
under grants FPA2004-00973 and FPA2007-62833, NSF (USA), US Department of 
Energy, NASA under the grant number NAG5-10842 and the helpful discussions within 
the network on direct dark matter detection of the ILIAS integrating activity (Contract 
number: RII3-CT-2003-506222). 

References 

[1] Peccei R D and Quinn H R, CP conservation in the presence of pseudoparticles, 1977 Phys. Rev. 
Lett. 38 1440 

[2] Peccci R D and Quinn H R, Constraints imposed by CP conservation in the presence of 

pseudoparticles, 1977 Phys. Rev. D 16 1791 
[3] Weinberg S, A new hght boson?, 1978 Phys. Rev. Lett. 40 223 

[4] Wilczek F. Problem of strong P and T invariance in the presence of instantons, 1978 Phys. Rev. 
Lett. 40 279 

[5] Kim J E, Weak interaction singlet and strong CP invariance, 1979 Phys. Rev. Lett. 43 103 

[6] Shifman M A, Vainshtein A I and Zakharov V 1, Can confinement ensure natural CP invariance 

of strong interactions?, 1980 Nucl. Phys. B 166 493 
[7] Dine M, Fischler W and Srednicki M, A simple solution to the strong CP problem with a harmless 

axion, 1981 Phys. Lett. B 104 199 
[8] Zhitnitskii A R, 1980 Yad. Fiz. 31 497 

Zhitnitskii A R, On possible suppression of the axion hadron interactions, 1980 Sov. J. Nucl. 
Phys. 31 260 (translation) 
[9] Amsler C et al (Particle Data Group), Review of particle physics, 2008 Phys. Lett. B 667 1 
[10] Duffy L D and van Bibber K, Axions as dark matter particles, 2009 [arXiv:0904.3346vl] 
[11] Zioutas K, Tsagri M, Semertzidis Y, Papaevangelou T, Dafni T and Anastassopoulos V, Axion 
searches with helioscopes and astrophysical signatures for axion(-like) particles, 2009 New J. 
Phys. (in press) [arXiv:0903.1807v4] 
[12] Sikivie P, Experimental tests of the "invisible" axion, 1983 Phys. Rev. Lett. 51 1415 

Sikivie P, 1984 Phys. Rev. Lett. 52 695 (erratum) 
[13] van Bibber K, Mclntyre P M, Morris D E and Raffelt G G, Design for a practical laboratory 

detector for solar axions, 1989 Phys. Rev. D 39 2089 
[14] Lazarus D M, Smith G C, Cameron R, Melissinos A C, Ruoso G, Semertzidis Y K and Nezrick F A, 

Search for solar axions, 1992 Phys. Rev. Lett. 69 2333 
[15] Moriyama S, Minowa M, Namba T, Inoue Y, Takasu Y and Yamamoto A, Direct search for 
solar axions by using strong magnetic field and x-ray detectors, 1998 Phys. Lett. B 434 147 
[hep-ex/9805026] 

[16] Inoue Y, Namba T, Moriyama S, Minowa M, Takasu Y, Horiuchi T and Yamamoto A, Search 

for sub-electronvolt solar axions using coherent conversion of axions into photons in magnetic 

field and gas hehum, 2002 Phys. Lett. B 536 18 [astro-ph/0204388] 
[17] Inoue Y, Akimoto Y, Ohta R, Mizumoto T, Yamamoto A and Minowa M, Search for solar axions 

with mass around 1 eV using coherent conversion of axions into photons, 2008 Phys. Lett. B 

668 93 [arXiv:0806.2230] 

[18] Zioutas K et al, A decommissioned LHC model magnet as an axion telescope, 1999 Nucl. Instrum. 
Methods Phys. Res. A 425 480 [astro-ph/9801176] 



Search for 14-4 keV solar axions with CAST 



22 



[19] Zioutas K et al (CAST Collaboration), First results from the CERN Axion Solar Telescope, 2005 

Phys. Rev. Lett. 94 121301 [hep-cx/0411033] 
[20] Andriamonje S et al (CAST Collaboration), An improved limit on the axion-photon coupling 

from the CAST experiment, 2007 JCAP 04 010 [hep-ex/0702006] 
[21] Arik E et al (CAST Collaboration), Probing eV-scale axions with CAST, 2009 JCAP 02 008 

[arXiv:0810.4482] 

[22] Paschos E A and Zioutas K, A proposal for solar axion detection via Bragg scattering, 1994 Phys. 
Lett. B 323 367 

[23] Avignone III ¥ T et al (SOLAX Collaboration), Experimental search for solar axions via coherent 
Primakoff conversion in a germanium spectrometer, 1998 Phys. Rev. Lett. 81 5068 [astro- 
ph/9708008] 

[24] Morales A et al (COSME Collaboration), Particle dark matter and solar axion searches with a 
small germanium detector at the Canfranc underground laboratory, 2002 Astropart. Phys. 16 
325 [hep-ex/0101037] 

[25] Bernabei R et al, Search for solar axions by Primakoff effect in Nal crystals, 2001 Phys. Lett. B 
515 6 

[26] Ahmed Z et al (CDMS Collaboration), Search for axions with the CDMS experiment, 2009 
[arXiv:0902.4693] 

[27] Krcmar M, Krecak Z, Stipcevic M, Ljubicic A and Bradley D A, Search for solar axions using 

5'^Fe, 1998 Phys. Lett. B 442 38 [nucl-ex/9801005] 
[28] Krcmar M, Krecak Z, Ljubicic A, Stipcevic M and Bradley D A, Search for solar axions using 

^Li, 2001 Phys. Rev. D 64 115016 [hep-ex/0104035] 
[29] Jakovcic K, Krecak Z, Krcmar M and Ljubicic A, A search for solar hadronic axions using ®^Kr, 

2004 Radiat. Phys. Chem. 71 793 [nucl-ex/0402016] 
[30] Derbin A V, Egorov A I, Mitropolsky I A and Muratova V N, Search for solar axions emitted in 

an Ml transition in ^Li* nuclei, 2005 Ptsma Zh. Eksp. Tear. Fiz. 81 453 [2005 JETP Lett. 81 

365 (translation)] 

[31] Namba T, Results of a search for monochromatic solar axions using ^^Fe, 2007 Phys. Lett. B 645 
398 

[32] Derbin AY et al. Search for resonant absorption of solar axions emitted in an Ml transition in 
5'^Fe nuclei, 2007 Pisma Zh. Eksp. Tear. Fiz. 85 15 [2007 JETP Lett. 85 12 (translation)] 

[33] Belli P et al, ''Li solar axions: Preliminary results and feasibility studies, 2008 Nucl. Phys. A 806 
388 

[34] Ljubicic A, Kekez D, Krecak Z and Ljubicic T, Search for hadronic axions using axioelectric 

effect, 2004 Phys. Lett. B 599 143 [hep-ex/0403045] 
[35] Bellini G et al (Borexino Collaboration), Search for solar axions emitted in the Ml transition of 

^Li* with Borexino CTF, 2008 Eur. Phys. J. C 54 61 
[36] Haxton W C and Lee K Y, Red-giant evolution, metallicity, and new bounds on hadronic axions, 

1991 Phys. Rev. Lett. 66 2557 
[37] Kaplan D B, Opening the axion window, 1985 Nucl. Phys. B 260 215 
[38] Srednicki M, Axion couplings to matter, 1985 Nucl. Phys. B 260 689 

[39] Mateu V and Pich A, V^^s determination from hyperon semileptonic decays, 2005 JHEP 10 041 
[40] Avignone III F T al, Search for axions from the 1115-keV transition of ^^Cu, 1988 Phys. Rev. 
D 37 618 

[41] Mallot G K, The spin structure of the nucleon, 2000 Int. J. Mod. Phys. A 15S1 521 [hep- 
cx/9912040] 

[42] Moriyama S, A proposal to search for a monochromatic component of solar axions using ^^Fe, 

1995 Phys. Rev. Lett. 75 3222 [hep-ph/9504318] 
[43] Turck-Chieze S et al, The solar interior, 1993 Phys. Rept. 230 57 

[44] Bahcall J N al. Solar models: current epoch and time dependences, neutrinos, and 
helioseismological properties, 2001 Astrophys. J. 555 990 



Search for 14-4 keV solar axions with CAST 



23 



[45] Raffelt G G, Astrophysical axion bounds diminished by screening effects, 1986 Phys. Rev. D 33 
897 

[46] Gondolo P and Raffelt G G, Solar neutrino limit on axions and keV-mass bosons, 2009 Phys. 

Rev. D 79 107301 [arXiv:0807.2926] 
[47] Kuster M et al, The X-ray telescope of CAST, 2007 New J. Phys. 9 169 [physics/0702188] 
[48] Abbon P et al, The Micromegas detector of the CAST experiment, 2007 New J. Phys. 9 170 

[physics/0702190] 

[49] Auticro D et al, The CAST Time Projection Chamber, 2007 New J. Phys. 9 171 [physics/0702189] 
[50] Luzon G et al. Background studies and shielding effects for the TPC detector of the CAST 

experiment, 2007 New J. Phys. 9 208 [arXiv:0706.1636] 
[51] Andriamonje S et al (CAST Collaboration), Search for solar axion emission from ^Li and 

D(p,7)3He nuclear decays with the CAST 7-ray calorimeter, 2009 [arXiv:0904.2103v2] 
[52] Chang HM et al (TEXONO Collaboration), Search for axions from the Kuo-Shcng nuclear power 

reactor with a high-purity germanium detector, 2007 Phys. Rev. D 75 052004 
[53] Derbin A V e< al. Search for solar axions produced by Primakoff conversion using resonant 

absorption by ^^^Tm nuclei, 2009 [arXiv:0904.3443vl] 



